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Abstract. The question is answered whether a 2 -shell-dynamos are able to produce a cyclic activity or not. Only kinematic 
dynamos are considered and only the solutions with the lowest dynamo number are studied without restrictions about the axial 
symmetry of the solution. The Q-effect is allowed to be latitudinally inhomogeneous and/or anisotropic, but it is assumed as 
radially uniform in the turbulent shell. 

For a symmetric a-tensor we only find oscillatory solutions if three conditions are simultaneously fulfilled: i) the a zz vanishes 
or is of the opposite sign as a^, ii) the a-effect is strongly concentrated to the equatorial region (i.e. it vanishes at the poles) 
and iii) the a-effect is concentrated to a rather thin outer shell. In the other cases almost always the nonaxisymmetric field 
mode S 1 possesses the lowest dynamo number which slowly drifts along the azimuthal direction. Also uniform but anisotropic 
Q-effect (a zz — 0) leads to the nonaxisymmetric solutions as it is confirmed by the Karlsruhe dynamo experiment. 
However, one of the antisymmetric parts of the a-tensor (not the vertical magnetic pumping) basically plays the role of a 
differential rotation in the induction equation. Using for the radial profile of this effect the results of a numerical simulation for 
the a-tensor of the solar convection zone (Ossendrijver et al. 2002), one indeed finds the possibility of oscillating a -dynamos 
even without the existence of real nonuniform plasma rotation so that they really can be called as pseudo-ai?-dynamos. The 
resulting butterfly diagram, however, proves to be of the antisolar type. The radial gradient of this pseudo-differential rotation 
determines the sign of the phase relation of B r and B<j, (see Sect. ^). 
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1. Introduction 

There are two observations which in the past led to a increas- 
ing interest in the solutions of the relatively simple a 2 -dynamo. 
The first one is the cyclic orbital modulation of close binary 
systems such as reported by Hall (1990) and Lanza & Rodono 
(1999) for RSCVn stars and the second one is the flip-flop phe- 
nomenon as reported recently by Tuominen et al. (1999) and 
Korhonen et al. (2001) for FK Coma stars but also in the sin- 
gle young dwarf LQ Hya (K2V, P Iot — 1.6 days, see Rice & 
Strassmeier 1998). Together with the highly nonaxisymmetric 
field configurations for very young cool dwarf stars reported 
by Jardine et al. (2002) (see Fig. [j]) there is increasing evidence 
that the traditional stationary and axisymmetric dipole-solution 
which is known since decades does not form the final truth. 

Oscillating a 2 -dynamos have been found to occur in the 
special circumstance wherein the a-effect changes rapidly 
in boundary layers (Radler & Brauer 1987; Baryshnikova & 
Shukurov 1987). In such cases, the period of the a 2 -dynamo 
depends strongly on the location of the a-boundary layer and 
is typically an order of magnitude or more smaller than the 
magnetic diffusion time across the dynamo-generation shell. In 
general, oscillatory dynamo behavior is produced by a combi- 
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nation of both the a-effect and the i?-effect. Kinematic mod- 
els of the oscillatory solar dynamo are always of the aJ?-type. 
Here, also motivated by the publication of Schubert & Zhang 
(2000), we discuss the question whether also a 2 -dynamos 
alone are able to exhibit oscillatory solutions. Schubert & 
Zhang (2000) considered a 2 -dynamos with the (academic) as- 
sumptions of completely homogeneous and isotropic a-effect 
and found oscillatory solutions for thin outer shells. In Section 
[3^2] we shall rediscuss this constellation under the extra condi- 
tion that only the lowest eigenvalue provides a stable solution. 
Then all solutions prove to be non-oscillating. 



The present paper only concerns kinematic dynamos on the 
basis of a-tensors derived by numerical simulations of MHD- 
convection. Nonlinearities in the mean-field dynamo equations 
are not considered. They are highly complicated as the majority 
of the resulting magnetic field configurations is nonaxisymmet- 
ric (see Moss & Brandenburg 1995; Kiiker & Riidiger 1999). 



Within the considered turbulent shell the a-effect is as- 
sumed as uniform in radius, there is no change of its sign. 
Stefani & Gerbeth (2002) have shown that magnetic oscilla- 
tions can occur if two shells with different signs of the a- 
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Fig. 1. Pole-on representation of the magnetic field of AB Dor. 
Courtesy M. Jardine. 

effect existj^j We do not follow this possibility in the present 
paper, here the a-effect has been considered as uniform in ra- 
dius within the turbulent shell. 



2. Basic equations and the model 

The model consists of a turbulent fluid in a spherical shell of 
inner radius ri„ and outer radius r out . A magnetic field is gen- 
erated in the shell by the a-effect (Steenbeck & Rrause 1966; 
Roberts 1972). In the shell, the turbulent magnetic diffusivity 
rjo is constant. For r > r out , we assume that there is a conduc- 
tor with large magnetic diffusivity rj out ; for r < r in , we assume 
that there is a conductor with high electrical conductivity, i.e. 
small magnetic diffusivity r) m . The induction equation is 



dB 

~dt 



curl(a o B — rj^curlB), 



(1) 



where B is the magnetic field, 77T is the turbulent magnetic 
diffusivity and a is the a-tensor. In no case the a-tensor has 
a too simple structure. For astrophysical applications it repre- 
sents the interaction of an anisotropic turbulence with a global 
rotation and a uniform magnetic field. There are also attempts 
to include the influence of a nonuniform rotation but this is 
beyond the scope of the present study. Here we are only inter- 
ested in the structure of the solution of a pure a 2 -dynamo, i.e. 
the rotation is assumed to be uniform. 

In Riidiger & Kitchatinov (1993) one finds the overall 
structure of the a-tensor as 

a im = - ai {G° i2°)5 vm - a 2 (G°/2,° n + fl?G° n )+ 
+a 3 (G° m f]f - G?<) - a 4 (G /2 )tf°< - je imk G° k . (2) 

1 The oscillating solution, however, seems to disappear for geome- 
tries where one sign of the a-effect dominates the other as it is the 



Here the unit vector f}° denotes the direction of the axis of 
the global rotation of the turbulence and the radial unit vector 
G° denotes its anisotropy. In almost all papers about a-effect 
dynamos the expression (Q) is reduced to its first term of the 
tensorial expression. We shall demonstrate in the present paper 
that only the inclusion of the remaining parts of the a-tensor 
reveals the variety of the solutions of the a 2 -dynamo and also 
solves the problem whether a 2 -dynamos can oscillate or not. In 
the whole paper the influence of the large-scale flow pattern is 
ignored. Then the remaining dimensionless number may be the 
dynamo number, C a = |ai|i?/^T with R as the stellar radius. 

It is interesting to consider the antisymmetric parts in 
the tensor (Q). A comparison with the induction term Si — 
tikmUkBm reveals that 7^ in the last term of (Jfy plays the role 
of a radial advection ("pumping") of the magnetic field. On the 
other hand, if formally a basic rotation with u = CI x x is used 
for the velocity field then one obtains £^ = (Q m Xi~- HiX^Bm. 
By comparison with (||) it follows that the -term in (g) ex- 
actly plays the role of a global (differential) rotation. More 
exactly speaking, in cylindrical coordinates (s, <fi, z) we find 
(a sz — «zs)/2s playing the role of an angular velocity, the gra- 
dient of which induces magnetic fields. In Section ^| we shall 
present the influence of the antisymmetric term 013 in (Q) for 
the question of the existence of oscillating a 2 -dynamos. 

3. The symmetric a-tensor 

We start in cylindrical coordinates with the symmetric part of 
the a-tensor, i.e. 

— a2sin0 

-aicosfl 1 .(3) 

— (ai + 2a 2 + a 4 ) cost 




As we shall demonstrate, the ratio 

ai + 2a2 + a 4 
a z = 



(4) 



will be of particular relevance for the resulting solutions. In 
spherical coordinates (r, 6, (j>) we have the general structure 



a. rr a r e 
n ~- 1 ae r aee 

a^tf, 

with a rr = — (ai + 2a 2 + a 4 cos 2 8) cos 9, a r g = a e 
a 4 cos 2 9) sin 9, agg = — (a x + a 4 sin 2 #)cos# and a 
— a 4 cos 9 so that for the ratio M) the expression 



(5) 

(a 2 + 

) 

(6) 



case in Fig. |2| 



a rr (pole) 
a^(pole) 

results where any of the two poles can be taken. 

3. 1. The a-tensor elements 



Ossendrijver et al. (2001, 2002) presented simulations for all 
the components of the a-tensor in spherical coordinates. Some 
of the results, which are relevant for the present discussion, are 
reported here. The simulations were done in a Cartesian box 
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Fig. 2. Numerical results for the Cartesian a-tensor compo- 
nents ctyy = a^cf, and a zz — a rr , measured in units of 
0.01 \fdg, as a function of depth in units of d from Ossendrijver 
et al. (2002) for the case of a box located at the south pole (run 
A00). The simulation domain consists of a thin cooling layer 
(z < 0), a convectively unstable layer (0 < z < d = 1), and 
a stably stratified layer with overshooting convection (z > d). 
The Coriolis number of the run is about 2.4, which is in the ap- 
propriate range for the bottom of the solar convection zone. For 
the other parameters and for a detailed description of the model 
we refer to Ossendrijver et al. (2001, 2002). The black curves 
are spatial and temporal averages; the shaded areas provide an 
error indication. Note, that the does not vanish at the pole. 

which is meant to represent a section from the lower part of 
a stellar convection zone, including a convectively stable layer 
underneath it. The Cartesian coordinate frame of the box corre- 
sponds to the spherical coordinates introduced above, such that 
the x-direction corresponds to the negative ^-direction, the y- 
direction to the ^-direction, and the z-direction (depth) to the 
negative r-direction. All simulations were done at the south- 
ern hemisphere of the star, and the angle between the vertical 
(radial) direction and the axis of rotation was varied between 
0° (south pole) and 90° (equator). Provided the rotation rate 
is sufficient to yield an a-effect, the depth dependence of a$g 
and has a typical shape, namely a negative sign in the bulk 
of the thick convection zone, and a positive sign in the thin 
overshooting layer. These features are expected for the south- 
ern hemisphere. The amplitudes of age and a^ increase with 
increasing angular distance from the equator up to a point close 
to the south pole, more or less consistent with the commonly 
assumed cos (^-function. Hence, the a-effect does not vanish at 
the poles. For weak rotation, the component a rr has a larger 
amplitude than the other two diagonal components, and it has 
the opposite sign. If rotation increases beyond a certain point, 
a rr as a function of depth developes multiple sign changes, un- 
like agg and O!00, and its amplitude falls behind that of the lat- 
ter. This is the rotational quenching of the vertical alpha effect 
reported in Ossendrijver et al. (2001). The symmetric compo- 
nent af s = a s zx is positive in the bulk of the unstable layer, 
and is generally larger in magnitude than a^ e ; it vanishes at the 
poles. 
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Fig. 3. Homogeneous and isotropic a-effect: The magnetic- 
field modes with the lowest value for C a for various shell 
thickness. All the solutions are non-oscillating. The curve is 
marked with the latitudinally mode number n for which the 
dynamo number is minimum. The oscillating modes found by 
Schubert & Zhang (2000) are located above the line with the 
lowest eigenvalues. 

3.2. Results for homogeneous a-effect 

There are references in the literature where the basic antisym- 
metry of the a-effect with respect to the equatorial midplane 
has been neglected (Radler & Brauer 1987; Schubert & Zhang 
2000; Stefani & Gerbeth 2000; Radler et al. 2002). The only re- 
alization of such models can only be imagined in technical ex- 
periments. Indeed, the paper by Radler et al. (2002) concerns 
the Karlsruhe dynamo experiment with a fixed helicity and a 
uniform flow field in vertical (z-)direction, so that a zz = 
is obvious. In order to demonstrate the differences between 
isotropic and anisotropic a-tensors also in the case of homo- 
geneous a-effect (i.e. cos 9 ignored), the following calculations 
are presented for the cases i) a zz = and ii) a zz — 0. While 
the first case seems to be an academic problem, the second case 
fits the situation in the dynamo experiment mentioned. 

We start with spherical models embedded in vacuum^ with 
outer a-effect zone located between x = x ln and x = 1. In this 
zone a-effect and eddy diffusivity tjt are assumed as radially 
uniform, below the convection zone there is a steep transition 
zone to the perfect-conducting interior (factor 100). Always 
the lowest dynamo numbers C a — aiR/i^ for the modes are 
given in the present paper. Any mode has an own dynamo num- 
ber, the mode with the lowest dynamo number is the preferred 
mode with the highest stability (Rrause & Meinel 1988). 

The problem of both homogeneous and isotropic a-effect 
(i.e. ai m — &o5i m with ao = const.) has been considered 
by Schubert & Zhang (2000). Our results do not comply with 
theirs if only the solutions with the lowest dynamo numbers are 
considered. In Fig. || the resulting minimum dynamo numbers 
C a and the associated latitudinal mode number n are given for 

2 Test computations for models embedded in a perfect conductor 
did not reveal basic modifications of our results in contrast to the re- 
mark by Brandenburg (1994) 
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Table 1. Dynamo numbers for anisotropic (a zz = 0) but uni- 
form a-effect (a = const., cf. Karlsruhe dynamo experiment). 
The boldface numbers mark the magnetic mode with the lowest 
eigenvalue 





AO 


SO 


Al 


SI 


0.25 


9.06 


8.40 


6.72 


6.72 


0.50 


10.84 (osc) 


10.64 (osc) 


9.50 


9.50 


0.75 


17.29 (osc) 


17.32 (osc) 


16.23 


16.23 


Table 2. The same as in Table |T] but for a ~ cos t 


^in 


AO 


SO 


Al 


SI 


0.25 


14.99 (osc) 


14.94 (osc) 


10.64 


9.97 


0.50 


15.7 (osc) 


15.5 (osc) 


11.8 


11.7 


0.75 


21.58 (osc) 


21.58 (osc) 


18.43 


18.43 



various inner shell radius x- ln . In opposition to the case of in- 
homogeneous a-effect (~ cos 9) the solution of the induction 
equation ([[]) are modes with a single latitudinal mode num- 
ber n. All the solutions can thus be characterized by the mode 
number n; the solutions with the mode number n and with the 
lowest eigenvalue C a are given in Fig. || - they all are non- 
oscillating. We do not find any oscillating a 2 -dynamo unless 
its a-amplitude was not the lowest one. As we take from the 
Fig. ||, however, it is clear that it is not enough to consider only 
the modes with n = 1 or n = 2. The thinner the (outer) a-shell 
the higher is the latitudinal mode number n of the modes ex- 
cited with the lowest dynamo number C a . Oscillating modes 
with low n are also existing but for thin shells their C a is never 
the lowest one.[] 

More interesting are the solutions with homogeneous but 
anisotropic a-tensor. The a-tensor has no zz-component, i.e. 
o-zz = in cylindric coordinates. Now the rotation axis is 
clearly defined so that it makes sense to ask for the axisym- 
metry of the solutions. Table [j] gives the results for uniform 
a-effect and Table || gives the results if the a-effect is anti- 
symmetric with respect to the equator, i.e. a ~ cos#. It was 
important to include the nonaxisymmetric modes into the con- 
sideration as they indeed possess the lowest dynamo numbers. 
Always the nonaxisymmetric modes with m = 1 dominate 
for all the models considered; and there are again no oscilla- 
tions. The axisymmetric solutions (mostly oscillating) which 
have been found by Busse & Miin (1979), Weisshaar (1982) 
and Olson & Hagee (1990) are probably not stable as the solu- 
tion with the lowest dynamo number are nonaxisymmetric and 
azimuthally drifting rather than axisymmetric and oscillating. 

There are no basic differences for the cases given in Tables 
[j] and |[ i.e. for a ~ const, and a ~ cosO (a zz = in both 
cases). The dominance of the m = 1 modes for anisotropic a- 
effect is a well-known result which has already been presented 
by Riidiger (1980) and by Riidiger & Elstner (1994). Here we 

3 The same remark concerns the analysis of Radler & Brauer 
(1987). 



have added its relevance also for the more simple case of ho- 
mogeneous a-effect. In the light of these calculations it is thus 
no surprise that the Karlsruhe dynamo experiment indeed pro- 
vides the nonaxisymmetric modes with m = 1 (Stieglitz & 
Miiller (2001). 

3.3. Results for inhomogeneous a-effect 

In the following we shall develop further the models with the 
equatorial antisymmetry of the a-effect as it results as the con- 
sequence of a global rotation of the considered spherical object. 
But we shall consider different a-profiles in latitude in order to 
simulate a possible concentration of the a-effect to the equator. 
We fix the latitudinal profile such as 

a oc sin 2A 8 cos 0, (7) 

where A is a free parameter describing the latitudinal profile of 
the a-effect. For A > the a-effect at the poles vanishes. It is 
more and more concentrated at lower latitudes for increasing 
value of A. Note, however, that the box simulations did not re- 
veal basic deviations of the latitudinal a-profile from the cos 6- 
law. In particular, at the poles the effect did not vanish. Insofar, 
the dynamo models with A > seem to be only of academic 
interest, but it is interesting to know that only for such models 
oscillating solutions appear with the lowest eigenvalues. 

The results for kinematic dynamo models with the (real- 
istic) equatorial antisymmetry are summarized in the Tables 
||. . .|] presenting the critical eigenvalues C a for different ex- 
tensor models and convection zones with various depths. Our 
notation is the standard one, i.e. Am denotes a solution with an- 
tisymmetry with respect to the equator and with the azimuthal 
quantum number m (see Krause & Radler 1980). The oscil- 
lating solar magnetic field mode mode is antisymmetric with 
respect to the equator, it is of AO type. 

Table |3| gives the results for isotropic a-effect (a z = 1). 
Again the boldface numbers represent the absolutely lowest 
eigenvalues C a indicating maximal stability. For the standard 
case with A = Table [3] provides the axisymmetric dipole AO 
as the stable mode. This result, however, strongly depends on 
the latitudinal profile of the a-effect. Already for A = 1 the 
preferred mode is nonaxisymmetric, i.e. Al and/or SI and for 
A = 2 and 3 we always find the Al mode as the preferred 
one. Our result is that except for the simplest latitudinal profile 
of the a-effect (i.e. the cos ^-dependence) the solutions are no 
longer axisymmetric, so that the work with axisymmetric codes 
only has a very restricted meaning. Although oscillating modes 
also appear they never have the lowest eigenvalue. Note that 
from all our latitudinal a-profiles only the cos 6>-dependence 
leads to finite a-values at the poles. 

For anisotropic a-tensor we even have a more clear sit- 
uation. Simplifying we worked with a zz = 0. The Table ^ 
presents the results. Always the solutions are nonaxisymmetric 
as always the oscillating axisymmetric modes possess higher 
eigenvalues. They occur, however, for the same a-anisotropy 
but for the thin convection zones with x- m = 0.8 and A > 
(Table |5|). For such a model where the a-effect does not ex- 
ist in the polar regions we find that the mode with the lowest 
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Table 3. Dynamo numbers C a for isotropic a-effect (a z = 1). 
The bottom of the convection zone is at x- m = 0.5 



A 


AO 


SO 


Al 


SI 







9.41 


9.42 


9.75 




9.76 




1 


28.8 (osc) 


28.7 (osc) 


26.7 




26.7 




2 


41.1 (osc) 


41.2 (osc) 


38.8 




39.2 




3 


51.9 (osc) 


52.0 (osc) 


49.8 




50.3 




Table 4. Eigenvalues C a for 




and x 


A 


AO 


SO 


A 1 


SI 





15.7 (osc) 


15.5 (osc) 


11.8 




11.7 




1 


35.0 (osc) 


33.8 (osc) 


32.7 




31.3 




2 


51.7(osc) 


49.1 (osc) 


49.3 




47.0 




3 


66.6 (osc) 


63.3 (osc) 


64.3 




61.4 





Table 5. Eigenvalues C a for a zz — and x- m = 0.8 



A 


AO 


SO 


Al 


SI 





26.3 (osc) 


26.3 (osc) 


23.2 


23.2 


1 


63.0 (osc) 


62.9 (osc) 


63.0 


62.9 


2 


89.8(osc) 


89.4 (osc) 


90.2 


89.9 


3 


112.9 (osc) 


111.8 (osc) 


113.4 


112.7 




Fig. 4. Pumping velocity 70 = j y , measured in units of 
0.01\/dg, as a function of depth in units of d from Ossendrijver 
et al. (2002) for the case of a box located at the equator (their 
Fig. 4). The tensor quantity 70 equals rJ?x- 

however, it transforms poloidal magnetic fields to toroidal mag- 
netic fields with a phase relation depending on the sign of 

In spherical coordinates the antisymmetric parts of the ex- 
tensor (^) can be written as 






sin 9 «3 





sin 9 «3 





-7 





7 






With the notation by Ossendrijver et al. (2002) one finds a 3 = 

-70/ sin6>, i.e. a 3 = -7 (equ). 



eigenvalue (i.e. the stable mode) forms oscillating axisymmet- 
ric magnetic fields (of quadrupolar equatorial symmetry, no 
dipoles). It should be underlined, however, that such a cyclic 
behavior seems to be a rather exceptional case as it only ap- 
pears if three conditions are fulfilled, i.e. 

i) the a-tensor must be highly anisotropic, 

ii) the a-effect must be concentrated to the equator, 

iii) the convection zone is rather thin. 

The latter condition remembers a similar finding of Radler 
& Brauer (1987) - but in contrast to their consideration this 
thin-shell condition is not a sufficient one for oscillating a 2 - 
dynamos. 

4. The antisymmetric a-tensor 

Let us now turn to the antisymmetric parts of the a-tensor. 
There is at first the a 3 -component in the tensor formulation 
(|2|). As mentioned above, the formation (ae r — a r g)/ (2r sin 9) 
formally acts as a (differential) rotation - so that in reality, if 
03 is not too small - all a 2 -dynamos can operate as (pseudo) 
a fl -dynamos which are known as oscillatory. We shall denote 
this virtual angular velocity by J?t with 

n T = -2*. (8) 

r 

Obviously, the ratio a^/ai will determine the ability of the 
a 2 -dynamo to operate as a (pseudo) aI?-dynamo. In any case, 



4. 1. Simulations 

From the simulations of Ossendrijver et al. (2002) the fol- 
lowing numerical results on 70 and j r can be mentioned. 
Longitudinal pumping (70) is generally the strongest pump- 
ing effect observed in the simulations unless the rotation axis 
is in the radial direction. It has a predominantly negative sign 
within the bulk of the unstable layer and the overshoot layer, 
which signifies that the mean field is advected in the retrograde 
direction. In most cases reported in Ossendrijver et al. (2002), 
there is also a thin layer near the top of the convection zone 
where the field is pumped in the prograde direction. The lon- 
gitudinal pumping effect is strongly dependent on latitude; it 
vanishes at the pole and peaks at the equator. 

The direction of vertical pumping Qy r ) is downward Qy r < 
0; i.e. 72 > 0) in the bulk of the unstable layer and in the 
overshoot layer. Near the top of the box there is a thin layer 
where the pumping is directed upwards (j r > 0). There is little 
dependence on latitude or on rotation. 

4.2. Results 

We start with the (unrealistic) case of radially uniform a 3 . Then 
after (||) is dfi/dr — a^/r 2 , so that for positive a3 there seems 
to be superrotation. Characteristic results are given in Table |[ 
For too small the nonaxisymmetric solution (Al) of the 
a 2 -dynamo is hardly influenced. For a3 ~ 10, however, we 

4 sign (BrB^) = sign (df2 T /dr) 
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a 




- 1 .5 t ^ 

0.70 0.80 0.90 1.00 

RADIUS 

Fig. 5. The same as in Fig. |2| but in the representation = 
T (r) (seeEq.g). 

Table 6. Dynamo numbers for anisotropic a (a z = 0) and 
uniform «3-effect. x; n = 0.5, «2 = 





AO 


SO 


Al 


SI 


3 


14.9 (osc) 


14.4 (osc) 


12.1 


12.0 


5 


14.3 (osc) 


13.7 (osc) 


12.9 


12.9 


10 


12.9 (osc) 


12.0 (osc) 


15.0 


15.1 



Table 7. Dynamo numbers for anisotropic a (a z = 0) and 
nonuniform a^-effect (see Figs. ||, ||). x m = 0.5, a-i = 



&3 


AO 


SO 


Al 


SI 


1 


15.1 (osc) 


14.8 (osc) 


12.1 


11.9 


2 


14.9 (osc) 


14.5 (osc) 


12.9 


12.8 


3 


15.0 (osc) 


14.6(osc) 


17.4 


17.4 


-10 


14.7 (osc) 


14.7 (osc) 


16.6 


16.6 



already find an oscillating quadrupole. Between a 3 = 5 and 
a?3 = 10 there is the transition of the kinematic dynamo from 
nonaxisymmetric drifting modes to axisymmetric oscillating 
modes. 

Obviously, the considered values of for the transition to 
a pseudo aJ?-dynamo are rather high. It is the non-uniformity 
of the J?t -effect which appears in the induction equation and 
may only allow the oscillations. We have thus applied to the 
model the radial profile in Fig. [| which results from real simu- 
lations (Ossendrijver et al. 2002). The profile is multiplied with 
the amplitude which is varied in order to find the various dy- 
namo solutions. The results are given in Table ^ and they are 
a surprise. Indeed, the amplification of the a 3 -effect only by 
a factor of 3 leads to the appearance of oscillating solutions 
as the most stable one. Again it shows quadrupolar equatorial 
symmetry. Dipolar symmetry only appears for the amplifica- 
tion factor of —10, i.e. the (say) aJ?T-dynamo resulting by the 
simulations of Ossendrijver et al. may lead to oscillating solu- 
tions but with quadrupolar symmetry. 



5. Discussion 

We have shown that for spherical configurations with outer tur- 
bulent shells the basic solution for a 2 -dynamos is a nonaxisym- 
metric mode drifting in the azimuthal direction. Oscillating ax- 
isymmetric solutions are seldom exceptions for a-effects that 
are 

i) highly anisotropic, 

ii) strongly concentrated to the equatorial region, 

iii) restricted to thin outer shells. 

In all other cases the mode with the lowest eigenvalue (which 
are considered here as the only stable one) is nonaxisymmet- 
ric and almost always of quadrupolar symmetry. The azimuthal 
drift is always of the same order as the oscillation frequencies 
are. Note that the box simulations by Ossendrijver et al. (2001) 
did not lead to vanishing a-effect at the poles. 

The same is true for a 2 -dynamos under inclusion of those 
anisotropic parts of the a-tensor which can be combined to an- 
tisymmetric components, i.e. Gifi m — G m fii. This term cor- 
responds to a (pseudo-) angular velocity, see Eq. (||), which 
induces electrical fields in the same way as a real differen- 
tial rotation of the plasma would do. The box simulations of 
Ossendrijver et al. reveal this pseudo- J? (called here SIt) as 
increasing outwards (see Fig. |]) so that in any case, together 
with the positive (in the northern hemisphere), an antiso- 
lar butterfly diagram results if the solution was axisymmetric 
and oscillating. The necessary dipolar solution, however, only 
exists if the J?t -effect is artifically amplified by a factor of 10 
(see Table |). 

We conclude that shell-dynamos without differential rota- 
tion on the basis of MHD simulations of rotating stellar convec- 
tion always possess nonaxisymmetric magnetic field configura- 
tions such as recently found for AB Dor. Oscillating solutions 
for a 2 -dynamos are revealed as rather exceptions. Our conclu- 
sion is that a cyclic stellar activity can always be considered 
as a strong indication for the existence of internal differential 
rotation. 
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